Landau level spectroscopy of surface states in the topological insulator Bio.91Sbo.09 via 

magneto-optics 



A. A. SchafgansQ B. C. Chapler, K. W. Post, and D. N. Basov 
Department of Physics, University of California, San Diego, La Jolla, California 92093, USA 

A. A. Taskin and Yoichi Ando 

Institute of Scientific and Industrial Research, Osaka University, Osaka 567-0047, Japan 

Xiao-Liang Qi 

Department of Physics, Stanford University, Stanford, California 94305, USA 
(Dated: February 21, 2012) 

We have performed broad-band zero-field and magneto- infrared spectroscopy of the three dimen- 
sional topological insulator Bio.9iSbo.o9- The zero- field results allow us to measure the value of 
the direct band gap between the conducting La and valence Ls bands. Under applied field in the 
Faraday geometry (k \\ H \\ CI), we measured the presence of a multitude of Landau level (LL) 
transitions, all with frequency dependence uj oc y/ll. We discuss the ramification of this observation 
for the surface and bulk properties of topological insulators. 



Bii_a;Sba^ [1-4 was among the first materials predicted 
to be a three-dimensional (3D) topological insulator (TI) 
[5 - 9^ : a material with bulk insulating properties that sup- 
ports conducting two-dimensional (2D) surface states. 
The appearance of nontrivial topological order is inti- 
mately tied to the band inversion that takes place as 
Bi is alloyed with Sb. Among the theoretical predictions 
now verified by experiment [TOHTB was for a free electron 
gas composed of spin-polarized quasiparticles existing at 
the 2D surface of the bulk 3D material, called surface 
states (SSs), formed when linear Dirac bands cross the 
Fermi energy {Ef). When a magnetic field is applied to 
the system, the SSs are thought to become gapped and 
spin-polarization reduced due to the time-reversal break- 
ing field. In this work, best illustrated by the results in 
Fig. |3j we demonstrate that all of the LLs observed in 
far-infrared magneto-optics obey Dirac-like dispersion as 
a function of applied field and are most likely due to op- 
tical transitions between LLs formed from the Dirac-like 
SS bands. 

The sample in this study was a large (~ 1 cm^) single 
crystal of Bio.91Sbo.09 cut along the bisectrix [211] plane 
[14]. This is the same plane (perpendicular to the [111] 
plane) in which quantum oscillations due to a 2D Fermi 
surface (FS) were observed in magneto-transport [14HT61, 
implying the presence of topological SSs. We measured 
near-normal incidence reflectance in the far-infrared (30- 
700 cm~^) as a function of temperature and applied field. 
In order to determine the zero-field optical constants, we 
measured reflectance between 30-8000 cm~^ and variable 
angle spectroscopic ellipsometry between 4500 - 45000 
cm~^. We extracted the optical constants by performing 
a Kramers-Kronig constrained variational analysis using 
refFIT software, based on a multi-oscillator fit of the re- 
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flectivity data anchored by the dielectric function mea- 
sured through ellipsometry [l7l [18] . These results were 
quantitatively similar to a full Kramers-Kronig inversion 
of the reflectivity data, which we show in the following 
figures. 



I. ZERO-FIELD RESULTS 



The zero-field data are presented in Fig. [T] showing re- 
flectance at four temperatures from room temperature to 
lOK (top panel). Unlike Bi2Se3 [19l[20], the plasma edge 
in Bio.9iSbo.o9, formed by the free electron response and 
characterized by the plasma frequency cjp, demonstrates 
substantial temperature dependence. At low tempera- 
tures, the lineshape of the plasma edge contains structure 
that may indicate the presence of multiple plasma fre- 
quencies and is indicated with an arrow in Fig. [l] With 
an Sb concentration of x = 0.09, the alloy is a direct- 
gap {Eg) semiconductor between the bulk La conduction 
band and the bulk Lg valence band [3 . This is evident 
in the real part of the optical conductivity cji (Fig. [l] 
bottom panel), where for photon energies greater than 
^ Eg -\- 2Ef, the onset of bulk interband transitions is 
observed: Eg -^2Ef = 470 ±10 cm~^, (indicated with an 
arrow in Fig. [l] bottom panel). The upturn of the con- 
ductivity spectra at low frequencies is due to the metallic 
Drude response. As temperature is lowered, the Drude 
response is significantly reduced, indicating a vanishing 
bulk metallic Fermi surface {ujp{T) is plotted in the in- 
set). However, the remnant Drude response at the lowest 
temperature shows that the material retains bulk metal- 
lic properties. 
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FIG. 1: Top. Reflectance vs. wavenumber in the far-infrared 
as a function of temperature, showing the Drude plasma fre- 
quency. The screened plasma frequency cJp(lOi^) is shown 
with an arrow. Top, inset A cartoon showing the bulk 
L bands, the Fermi energy Ef defined from the bottom of 
the conduction band, the surface states formed when the sur- 
face Dirac bands cross Ef, and the necessary photon energy 
Eg + 2Ef to create an interband transition. As far as we 
are aware, there have been no direct measurements of the SS 
bands along the CI plane in Bii-^Sbcc. Transport measure- 
ments suggest there is at least one SS with Dirac dispersion, 
and in principle, there can be more (3, 5, ...). Bottom. The 
real part of the optical conductivity ai{uj) as a function of 
temperature, illustrating the decreasing Drude response at 
lowered temperature. The optical gap Eg + 2Ef is shown 
with an arrow, and the conductivity rapidly increases at fre- 
quencies above the optical gap. Bottom, inset. The Drude 
plasma frequency 0Jp{T) obtained from modeling ai{uj) with 
the Drude-Lorentz model. The strong temperature depen- 
dence of 0Jp{T) is likely the consequence of thermal activation. 



A. Location and size of Fermi Surface 

Taskin and Ando [21] measured the volume of the 
bulk FS, formed by a set of three ellipsoids located at 
the L-points of the Briullioin zone, to be n = 8.1 ± 
0.2 * lO^^cm"^. The measured Hall coefficient implies 
a similar electron concentration of n = 1.8 * 10^^ cm~^ 
[2T1 . These data were obtained with samples from the 
same rod as the crystal used in this study and we there- 
fore expect the properties to be quantitatively similar. 
The Drude model of the lowest temperature conductiv- 
ity should produce comparable values between optics and 
transport only if the Fermi energy is located at an equiv- 
alent location in the bulk conduction band. This is be- 
cause the plasma frequency of the Drude response is re- 
lated to the number density of free carriers as: uo'^ = 
4 i\ne^ Ivff . Using the plasma frequency as measured 
via optics, uj^ = 1485 cm~^, the resulting band mass is 
m* = 0.0073me. The plasma edge seen in optics is dom- 
inated by the shortest of the semimajor axes of the FS 
ellipoid, measured to be /c/ = 2.3 * 10^m~^. The rela- 
tionship Vf = hkf/m'^, where is the Fermi velocity, 
should provide the band velocity of the linear bulk con- 
duction band along the fastest direction and therefore, 
the Fermi velocity in the fast direction of the bulk FS is 
Vf =3.65*10^ m/s. 

Next, we invoke Ef = hkfVf in order to determine 
the separation between the bottom of the conduction 
band and the Fermi energy. This gives Ef = 5.37meV 
= 42.5 cm~^ above the bulk band gap, into the bulk 
conduction band. Based on the optical conductivity, 
since Eg + 2Ef ^ 470± 10 cm~^, the magnitude of the 
bulk band gap between the La and Lg bands must be 
Eg ^ 387± 10 cm~^. This is in excellent agreement 
with the data in Fig. 1 as well as previous photoemis- 
sion studies. {Eg is very sensitive to Sb content and 
above x=0.04. Eg grows with increasing Sb content [3]. 
In slightly higher doped Bio.gSbo.i, a lower bound of 
£;^=50meV=403cm-^ was found using ARPES [10 .) A 
cartoon of the bulk L bands and the surface states is 
plotted in Fig. [TJ inset, showing the required incident 
photon energy to induce an interband transition between 
the bulk La and Lg bands. Also shown is a possible 
schematic of the SS bands [lOl El [22] , illustrating that 
multiple SSs may still exist when the Fermi energy lies in 
the bulk conduction band. We explicitly note that most 
of the cited works have been performed or calculated for 
the [111] surface, whereas our present work is performed 
on the [2ll] surface. The surface state band structure 
along this surface has not presently been determined by 
photoemission studies. 



II. IN-FIELD REFLECTANCE DATA AND 
MODELING 



Turning to the in- field data. Fig. [2] shows reflectance in 
O.IT increments at lOK (top) and lOOK (bottom). The 
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FIG. 2: Reflectance vs. wavenumber and applied field H || 
CI II A;, at T = lOK (top) and lOOK (bottom). Each curve 
represents O.IT increment in applied field, going from 0-3T at 
lOK and 0-2T at lOOK. 



magnetic field was applied parallel to the bisectrix axis 
(CI), perpendicular to the bisectrix plane (C2-C3 plane), 
and parallel to the A:-vector of the incident electric field 
(Faraday geometry) in order to access the previously ob- 
served 2D FS deduced from quantum oscillation (QO) 
measurements [14]. The large LL absorptions become 
visible above ~ 0.3T, which likely correspond to bulk 
states, as well as other smaller field-dependent features. 
Models of the lOK data were used to extract values of 
the cyclotron resonance (CR) (see appendix). The mod- 
els do not allow for the observation of the much weaker 
LLs due to the SS bands. Therefore, extracting LLs from 
the raw reflectance via modeling is not possible for any- 
thing other than the most prominent LLs because the 
surface state charge density is very small and therefore 
produces exceedingly weak features in reflectance. 

Instead of modeling the data, differentiating the re- 
flectance with respect to the applied field is a far more 
sensitive technique and allows access the subtle features 
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FIG. 3: Contour plots showing the derivative of the re- 
flectance with respect to magnetic field, dR/dH, at T = lOK 
(top) and lOOK (bottom). LLs show up as peaks that disperse 
to higher wavenumber with applied field. 



corresponding to the 2D FS that are otherwise dominated 
by the bulk properties. Once we perform a derivative of 
the reflectance with respect to the applied field, dR/dH, 
the weakest LLs become quite evident. Indeed, these 
features are so small that they would have gone unob- 
served without utilizing this technique. Fig. [3] shows the 
dR/dH surface contour at at both T= lOK (top) and 
lOOK (bottom), from uj= 40-680 cm"^ and H= 0-3T (at 
lOK) and i^=0-2T (at lOOK) with O.IT increments. The 
value of dR/dH corresponds with the color scale. The 
intermediate field data are interpolated. As can be read- 
ily seen, there are multiple LLs, as well as features (dark 
red) that correspond to the field-dependent behavior of 
the bulk Drude plasma edge. All of the observed LLs 
display uj cx dependence. 



4 



III. EVOLUTION OF LANDAU LEVELS IN 
MAGNETIC FIELD 

A classical picture explaining the origin of the CR is 
that when a magnetic field is applied to a metal, free 
charge carriers become bound in orbits around the field 
lines, provided the mean-free path is sufficiently long. 
Such bound charge carriers create a resonance in the op- 
tical refiectance, dependent upon the strength of the ap- 
plied field, the carrier's effective mass and number den- 
sity. Since the band dispersion near Ef is quadratic in 
most materials, the frequency of the CR and resultant LL 
transitions are ujc oc H. As was noted early on, LLs in ele- 
mental bismuth have a non-linear field dependence due to 
the linear band dispersion of the Ls and La bands |2 
[25j. Recently, the bulk Lg valence band in Bio.gSbo.i has 
been demonstrated to be linear using ARPES [10] , in ad- 
dition to the observation of linear Dirac bands that form 
SSs. 

In order to describe the energy spectrum of the LLs 
due to linearly dispersive bands, one must use a model 
derived from Dirac theory [23, 26 . For the bulk bands, 
the allowed energy states disperse as: 



E = ±{ 



El 



eHh 



(1) 



where /i^ is the Bohr magneton, nic is the experimental 
cyclotron mass, n is the LL index of the bulk states, go 
is the experimental ^-factor, and the first ± selects a LL 
in the conduction (+) or valence (-) band while the ± on 
the last term selects the spin state. This is different from 
the SSs in a three dimension topological insulator [27] . 
where the expression for the surface state LLs is: 
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Here, gg is the effective magnetic factor of the surface 
electrons, m is a correction to the effective Drude mass, 
N is the LL index of the SSs, and ± selects a LL above 
(+) or below (-) the Dirac point [27 . At low fields, the 
LLs disperse ex ^/H, but as the field increases the linear 
terms will dominate. 

In Fig. [4] we plot the single particle LLs for both the 
bulk (top) and SSs (bottom). There are several impor- 
tant features to be learned from these figures, as far as al- 
lowed optical LL transitions. First, as can be seen in the 
top panel, any interband LL transitions of the bulk states 
will extrapolate to a value Eg at zero-field, while bulk 
intraband transitions will extrapolate to zero-frequency 
at zero-field. Second, bulk intraband transitions are not 
allowed once all of the bulk conduction LLs rise above 
Ef. (Using the parameters mentioned in the figure cap- 
tion, this takes place near 1.5T. Of course, using a larger 
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FIG. 4: Single particle LLs as a function of frequency for 
the bulk (top) and Dirac (bottom) bands, showing the n = 
0-2 bulk LLs and the n — 0-3 Dirac LLs. The green hori- 
zontal line illustrates the location of the Fermi energy with 
respect to the single particle LLs (^/= 45 cm~-^) as extracted 
from the zero-field data, as discussed above. The bulk LLs 
are spin- split, with the size of splitting determined primar- 
ily by the strength of spin-orbit coupling, which bears on the 
value of go. The red (blue) lines are spin up (down) LLs. 
Values used to determine the LLs are, for the bulk states: 
77ic=0.0073me, ^^=385 cm""*^, go=0.5/mc. For the SSs, we 
used: m=l, ^s=80, Vf =8.5:^.10^ m/s. 



(smaller) effective mass means bulk intraband transitions 
are allowed to higher (lower) fields.) Third, all of the LL 
transitions due to the SSs (Fig. |4j bottom) will extrapo- 
late to zero-frequency at zero-field. Fourth, the 0th LL is 
almost field independent for the parameters we use (see 
Fig. [4] caption), and there is no spin-splitting. Fifth, LL 
transitions will not become allowed until the LLs cross 
Ef^ and therefore the onset of optical transitions should 
be evident in the data. We note that the bulk band gap 
Eg and Fermi level Ef may not be field-independent, but 
any such field dependence should be a higher-order cor- 
rection and will not significantly impact the results we 
show here. Also, the placement of the Dirac point, from 
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dR/dH map for BiSb at 10K, H || C1 
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FIG. 5: The theoretical magnetic field dependence, plotted 
as a dR/dH contour, of a plasma edge of a metallic semi- 
conductor similar to what we observe in Bii_ccSbcc. The quan- 
titative similarity with the data is suggestive. 



where the surface state LLs disperse in field, may not be 
located in the center of the bulk band gap. 

Additionally, we show that a simple model of the mag- 
netic field dependence of the Drude plasma edge in metals 
[26] can reproduce the behavior of the most prominent 
field-dependent features. Figure [5] demonstrates how the 
cyclotron active and inactive modes in magnetic field lead 
to a splitting of the plasma edge, where the cyclotron ac- 
tive mode moves the plasma edge higher in energy while 
the inactive mode suppresses the plasma edge. Quali- 
tatively, such behavior is present in the data we show. 
However, more detailed analysis is needed in order to 
make definitive quantitative conclusions. 



IV. DISCUSSION 

Based on the information gained from the single par- 
ticle LLs, we can make several immediate conclusions 
about the optical LL transitions in Bio.91Sbo.09 presented 
m Fig. |3l First, none of the observed LLs extrapolate to a 
zero- field value equal to the bulk band gap Eg. Therefore, 
based on eq. [l] none of the LLs are due to bulk interband 
transitions between LLs in the La and Ls bands. Second, 
as illustrated in Fig. [4j if there were any intraband LL 
transitions due to the bulk conduction bands present, 
these transitions should disappear by ~ IT. Therefore, 
none of the LLs can be understood to be due to the bulk 
bands within the paradigm presented by eq. [l] Third, 
the six highest energy LLs are observed to "turn on" at 
finite field, consistent with the notion that once the LLs 
cross above Ef, transitions become allowed. 

A logical place to begin to quantitatively understand 
the LLs is to assume that the same surface state observed 
in quantum oscillation measurements on nearly identical 
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FIG. 6: Theoretical LL transitions of the form in eq|2] over- 
laid on the data at lOK. We have included three different 
Fermi velocities in order to explain all of the observed LLs. 



As labeled, 



: 8.5xl0^m/s (black symbols), v'j 



■ 6xl0^m/s 



(white symbols) and Vf = 5x10 m/s (red symbols) and the 
0-1 (dashed line), 1-2:2-1 (crosses), and 2-3:3-2 (open circles) 
transitions are given. We did not correct for the Drude ef- 
fective mass (m=l) and used a very moderate g-factor of the 
surface states (^s = 120). 



samples is also contained in our data. Results from QO 
measurements determined the presence of a 2D FS with 
Vf = 8.5xl0^m/s for very low fields (<1T) and at low 
temperature [14] . In Fig. [g] we show the predicted LL 
transitions due to a Dirac band with a Fermi velocity of 
vj^ = 8.5xl0^m/s for the n=0 to n=l transitions (blue 
lines), 1-2:2-1 (black lines), and 2-3:3-2 (red lines) tran- 
sitions. Therefore, we confirm previous measurements of 
the 2D FS. 

Additionally, we find that by applying eq. [2j the num- 
ber of LLs is greater than what can be explained with a 
single surface Dirac band. That is, in order to explain 
the data, we must include transitions that would exist in 
the presence of two additional Dirac bands with Fermi 
velocities determined by fitting the LLs: v'j = 6xl0^m/s 
and vj = 5xl0^m/s. This presents one possible physi- 
cal interpretation of the origin of the weakest LLs ob- 
served. In this scenario, it is interesting that the 0-1 
transitions from v'j and vj are either not present or com- 
pletely overwhelmed by the behavior of the bulk plasma 
edge. ARPES studies of the [211] surface are required 
to determine the number of Dirac bands, and should be 
able to use our LL observations as a guide. 

The frequency at which these transitions begin pro- 
vides a clue as to the separation between the upper 
(+) and lower (-) LLs and ultimately provides a mea- 
sure of the location of the Dirac point {Eoirac) defined 
with respect to the bottom of the conduction band: the 
optical LL transitions of the SS bands will appear at 
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uj=2{Ef—EDirac)i and will extrapolate to zero-frequency 
at zero-field. At r=10K, the higher order transitions be- 
gin between 425-500 cm~^ while the 0-1 LL transition is 
seen to begin at an energy equal to 1 /2 of the higher order 
transitions: {Ef — E]jirac)' For instance, at T = lOK, 
the 0-1 LL vj transition first appears near 225 cm~^. For 
Ef = 42.5 cm~^, and therefore Fjjjirac is located ^182.5 
cm~^ below the bottom of the bulk L^ band, ^11 cm~^ 
above the center of the bulk band gap. 

Interestingly, we do not observe bulk LL transitions 
across the gap [28 . This is somewhat surprising, given 
the sensitivity of our method for detecting LLs. One pos- 
sibility is that the bulk transitions have short lifetimes, 
much shorter than the lifetimes of the SSs. The width 
of the Drude peak is one clue to the expected bulk LL 
lifetime, and roughly has a value of F = 80-120 cm~^, 
depending on modeling parameters. On the other hand, 
the lifetime of the SSs as determined from the width of 
the LL resonances is F = 10-30 cm~^, implying the SS 
carriers have lifetimes as much as 10 times longer than 
the bulk carriers. Transport measurements suggest a 
similar picture, where the mean free path of the bulk 
carriers was 16 nm while the SS carriers was 150 nm 
[15]. Our data, therefore, is consistent with the notion of 
scattering-protected SSs, while the bulk states are more 
strongly scattering. 

In conclusion, we have detailed a magneto-infrared 
study of the topological insulator Bio.91Sbo.09, demon- 
strating the presence of a multitude of Landau levels with 
UJ (X \fR dependence [9 . After careful analysis, the LLs 
that we observe can only be understood to arise from 
Dirac-like bands within the bulk band gap. We there- 
fore provide direct evidence for the existence of topolog- 
ical surface states in Bio.91Sbo.09, coexisting with a bulk 
metallic Fermi surface. Future temperature and doping 
dependent LL studies as well as gated structures where 
the chemical potential can be tuned into various energy 
regimes of the band structure will help to further eluci- 
date the origins of the observed phenomena we present. 

We thank Gil Refael and Doron Bergman for enlight- 
ening discussions. Funding was provided by DARPA and 
FENA. 



V. APPENDIX 

A. Reflectance models 

In order to extract the cyclotron resonance due to LL 
transitions, we attempted to model the far-infrared re- 
flectance using refFIT [17]. The resultant models are 
shown overlaid on the corresponding field-dependent re- 
flectance data in Fig. [7| where we point out the value of 



the modeling result for the primary cyclotron resonance 
with an arrow. The models do not allow for the observa- 
tion the much smaller LLs due to the SS bands. At fields 
closer to 3T, fitting converged on the presence of three 
CR modes, however they did not display an intelligible 
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FIG. 7: Reflectance in O.IT increments, overlayed with the 
corresponding model (red line) from RefFIT. 



systematic behavior. This may be due to the difficulty 
of fitting the raw reflectance and as such, we have not 
drawn conclusions based on these models. In Fig. |8| we 
plot the values of the CR as determined by modeling the 
raw reflectance (circles), overlaid on the dR/dH contour. 
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